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and moreover it shows the physical meaning of the parameters 
used in the simulation. 
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1 Introduction 

In a recent paper [T] a lattice gauge theory for massive SU(2) Yang-Mills 
in dimension four has been proposed. The completely gauge independent 
analysis shows the presence of a transition line in the parameter space 
(to 2 ~ which indicates a phase transition of second order for large 

/3, starting from (3 ~ 2.2. The /3 — > oo limit value of /3m 2 is in good agree- 
ment with the critical point of the nonlinear sigma model in four dimensions 
quoted in Ref. [2], namely 0.65 ± 0.04 (see Fig. EJ). 



In this paper we investigate on the small coupling, i.e. large /3, behavior 
of some global observables such as the "energy of the vacuum" and the order 
parameter 

where ./V is the number of sites and Zl the partition function Q It should 
be noticed that in Ref. pQ the order parameter is defined as 

1_ ^2^' (2) 

in order to put in evidence its behavior for to 2 = and to 2 — > oo. Here 
these properties are not relevant. 

In this work we discuss two important questions: i) to understand the 
phenomenological meaning of the parameters /3 and to 2 (in the naive con- 
tinuum limit, a — > , /3 = 4<? -2 and too -1 is the mass of the vector mesons); 
ii) to compare the lattice regularization of the theory with the continuum 
Minkoswkian formulation of massive Yang-Mills theory proposed in Refs. [3] 
and [3]. Regarding i) we want to find the scales necessary in order to inter- 
pret the otherwise dimensionless parameter to. For the point ii) we assume 
that the finite lattice artifacts are not relevant quantitatively for the global 
quantities of eq. (P). 

The transition from euclidean to minkoswkian quantities is performed 
in the usual conventional way: exp^S^) is the weight in the path integral 



3 -m for Minkowski, -e for Euclidean and -l for Lattice. 
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and the arrow of time (i.e. anticlockwise Wick rotation) is chosen in order 
to match the edge-of-the-wedge theorem. 

In the subtraction procedure [5] - [7] a scale A of the radiative correction 
is introduced, then it is natural to look for a correspondence implemented 
by the mapping 

g 2 = f(/3,m 2 ) 

Ma = t(f3,m 2 ), (3) 

where the parameters g, M are the field theory coupling constant and mass 
respectively, while a is the "lattice spacing". Eqs. ([3]) show that the mapping 
of the continuum field theory onto the lattice implies a choice of a surface 
in the three-dimensional space spanned by the dimensionless parameters (g, 
Ma and Aa). The convenience of these choices might be sustained and 
tested on the physical observables to be compared. 



2 Massive Yang-Mills in Field Theory (Minkowski) 

In the continuum quantum field theory the classical action is 

S m = —p- Jd D x[- z G aiMU [A]G^ [A] + — (A afl - F a ,f) . (4) 

The path integral functional Zm is obtained by integration over the fields 
A^ and f2 £ SU(2), where (r are the Pauli matrices) 

= 0o + iTj4>j, = iUd^. (5) 

The constraint on 

I = + (6) 

renders the theory in eq. (|3D nonrenormalizable. The procedure of subtrac- 
tion of the infinities in perturbation theory is described in detail in Ref. [3]. 
Here we account only for the final practical rule: only pure pole subtraction 
is performed on the dimensionally regularized Feynman amplitudes normal- 
ized by 

A( 4 " D U. (7) 
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Figure 1: The Goldstone boson lines are dashed. The Faddeev-Popov prop- 
agators are dotted. 

No additional finite adjustment is allowed, e.g. on-shell normalization. 

We use Landau gauge because it is very useful in a massive theory, 
where the massless modes are all unphysical. The absence of the massless 
modes contributions provides a good check for the calculation of the physical 
observables. 

The "energy of the vacuum" per hypercube a 4 is then given by the path 
integral mean value ()m 

£ M (9, M, A, a) = ( - ^ ( - \G a , v [A]G^[A] + . . . ) ) ^ (8) 

where . . . stand for the counterterms and for the terms of the classical action 
depending on the unphysical fields as Q and Faddeev-Popov ghosts c, c. The 
full expression is given in Ref. [3]. 

The O(g ) of Em is zero. This result can be proved formally by applying 
the operator g-^ to Zm before and after the rescaling of the field A^ — > gA^. 
The 0(g 2 ) part amounts to the evaluation of the graphs in Fig. (TJ i.e. after 
rescaling 

£m = { a % lf q (- \Go^[A]G^[A] + ^Al . . . )) m . (9) 



2.1 Two Loop Amplitude 

A straightforward calculation gives the Feynman amplitude properly nor- 
malized [5] by the factor A^~ D "> 



1M 



' 4 2 

-a g 



3(D - ^)M 2 I 2 [M] - (D 2 -AD + ^Wm] 



4 / 



(10) 
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The integrals are denned by 



I 2 [M] = 
A [M] 



d D c 



d D p 



A (8- 



2D) 



d (27T)D J m (27T)D (g2 _ M 2 )(p 2 _ M 2 )[{q + p) 2 _ M , 

dPq A^ D ) 



(11) 



' M (^) D (q 2 -M 2 )' 

It should be noticed that all contributions by the massless modes cancel out 
(no contributions from the unphysical sector). 
The tadpole integral is 

M 2 



A [M] = i 



(4tt) 



D-4 
1 



In A 



(D-4) In A + l 



7T~ 

~6 



where 



Consequently one gets 



A 



4vrA 2 ' 



+ 0((D-4) 2 ), (12) 



(13) 



A [Mf 



M 4 



TF> + n jlnA + l + ^ r +21n 2 A 

(D - 4) 2 D - 4 6 



(14) 



The two- loop integral can be found in the literature (see Ref. [8]} 



h[M] = 3 



M 2 
(4vr) 



+ 



{D-4) 2 D-4 



1 - 2 In A 



7T 



lnA + — + -+ln 2 A +—CU- 



12 



V3 



7T 



(15) 



where the Euler-Mascheroni constant and the Clausen function are 

7 = 0.5772156649 



Ck 



1.014941606. 



(16) 



Finally the contributions of the graphs in Fig. [T] add up to 

69 1 1 



-1M 



3 4 2 M 4 
-2 a9 j4^ 



91 69 , . 

1 in A 

2 (L>-4) 2 (D-4)\ 4 2 



+ — In A - — In 2 A - — tt 2 + 



163 91 

T 



23 
16' 



2V3 V3 



(17) 
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2.2 Contribution of the Counterterms 

In order to evaluate the contribution of the counterterms we use the results 
of Ref. gj eq. (29) 



f« = J d D x 7 (1) (x) 
A (D-4) x 



d D x 



(4tt) 2 D-A 
Then the contribution of the counterterms is 
.a 4 g 2 M 2 D-169 



^ (d,A au d»A: - {dA a f) + \M 2 Al] .(18) 



&2M 



(4vr) 2 D-44 



■A [M 



3 a 4 g 2 M 4 
2 (4vr) 



69 



1 / 69, \ 23, A 

„ H 23 H In A H In A 

L>-4 2 D-4V 2 / 2 



69/, 2a , vr 2 

H f In A + H 

8 V 6 



(19) 



The final result for the energy of the vacuum to order 0(g 2 ) is 



4„2 



£m = £im + £im = — 2° 9 
69 



M 4 



69 



1 



+ 



183 1 



2 (L>-4) 2 4 (£>-4) 

137, . 47 99 /7T 

— In A H In A 1 -Ch - 

8 4 4 2^3 3 



(20) 



Notice that the singular part has no dependence on A. The finite part is 
then 



-M 



3 4 2 M 4 

-« 9 TT=U 



(4tt)^ 



69, 2 . 137 47 99 _ /vr 

— In A H In A 1 -Cl 2 - 

8 4 2^3 



• (21) 



2.3 The Order Parameter 

The order parameter at 0(g°) is given by 

M 2 , . „ , 9 \(0) _ 3a 4 M 4 



r(0) 
'M 



/ m 2(4vr) 2 \L>-4 



31nA + 2), (22 



where the path integral mean value is taken with the free part of the action. 
The 0(g 2 ) can be easily evaluated with a trick. We use the fields scaled in 
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such a way that the coupling constant appears only in the interacting part 
of the action. Then we have the identity 

A2) _ g d (2 ) _.M 2 d g 8 



■M 



2 dg ~ M " N dM 2 2 dg 
Thus we use the expression in eq. (f20|) 



In Z M 



-Af" 







dM 2 



£m- (23) 



r(2) 
~M 



3 2 a 4 M 4 
off 



(47T) 



69 



1 



+ 



183 1 



(L>-4) 2 2 D-4 
69 , o . 205 , . 43 



69 , 2 . 2U5 . 43 00 /-^, /7T 
— In 2 A + — In A + — + 33V3C/2 - 
4 4 4 \3 



(24) 



We notice that once again the singular terms depending on A cancel out in 
the final amplitudes. 



3 Euclidean-Minkowskian Correspondence 

The correspondence is established by requiring that the observables have 
the same value. The mapping is obtained by the substitution 

£4 = — ixo 

A A = iA (25) 

and by performing an anticlockwise Wick rotation on the xq integration (to 
match the statement of the edge-of-the- wedge theorem). The generating 
functionals are obtained by summing over the field configurations with the 
weights 

e -S E 

e iSM (26) 

for the Euclidean () e and Minkowskian Qm mean values, where Se is the 
Euclidean action obtained from Sm by using the mapping in eq. (|25|) . Then 
we use the correspondence 

([...Hx 3 )...A, k (x k )...)) M \ ((...ol. )...,l / , i (.r,)...)) / ,i2T) 



A = -iAi 



In the Minkowski case the field theory is made finite by the procedure 
briefly outlined at the beginning of Sec. [2j For the Euclidean formulation 
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we want to investigate the possibility to approximate the amplitudes by a 
gauge lattice model given by an appropriate action Sl introduced in Rcf. [1 
and discussed later on. In fact the paper is devoted to this comparison by 
considering the global quantities Sl and £l in lattice gauge theory in the 
limit of weak coupling (/? S> 1), where perturbation theory can be used for 
the theory in the continuum. 



4 The Lattice Model 

The action in eq.Q is invariant under <?l(x) G SU(2)l local-left and gn E 
SU(2)r global-right transformations. On the lattice variables we have 

[ n'(x) = g L (x)Sl(x) [ n'(x) = n(x)gl , no . 

SU(2) L l K ' 3 V ' V ' . , SU(2) R t y ' y jyR .(28) 

[ U'(x,fi) = g L (x)U(x,n)g L (x + fi) [ U'(x, fi) = U(x, n) 

The lattice model is constructed by assuming the same invariance properties. 
The nearest neighbor interaction is required to map naively into the action 
@ for zero lattice spacing. The link variable is taken to be 

U(x,n) ~ exp(-iaA^(x)). (29) 



Thus the action is 

S L =| D\e^2Tr{l - U u } + ^m 2 <Ke^Tr{l - n(x)^U(x, n)Q(x (30) 

□ X\i 



where the sum over the plaquette is the Wilson action [9] and the mass term 
has the (Euclidean) continuum limit 



2' 



xp, 

M 2 

9 2 



' ; M 2 a 2 9\e J] Tr { 1 ~ ^(x)^U(x, n)Sl{x + //) J 

Xfl 

J d 4 x Tr{(A^-ind^) 2 }. (31) 

5 The Lattice Simulation 

The partition function is obtained by summing over all configurations given 
by the link variables and the <SC/(2)-valued field £1 

Z L [f5,m 2 ,N] = ]T e~ s \ (32) 
{u,n} 
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Figure 2: The transition line. The arrow marks the position of the end point. 
In the figures of the present paper the statistical errors are not displyed since 
they are too small to be shown. 

where iV is the number of sites. 

In principle the integration over Q(x) is redundant, since by a change of 
variables (Un(x, fi) := fl(x)^U(x, /i)0(x + /i)) we can factor out the volume 
of the group. Zi\f3,m 2 ,N] becomes 

[E] ^^-^(9te^Tr{l-U n } + m 2 9te^2Tr{l-U(x,Li)}). (33) 
{a} {u} □ x ^ 

However in eq. (|32p we force the integration over the gauge orbit Uq by 
means of the explicit sum over Vt. In doing this we gain an interesting 
theoretical setup of the model; in practice, our formalism is fully gauge 
invariant. Moreover by forcing the integration over the gauge orbit Uci 
we get results which are less noisy than those obtained by using only the 
integration over the link variables in eq. (|33p . 

5.1 Numerical results 

The numerical analysis of the model pQ shows the existence of a line (Figj2]) 
where the functions 8 and (£ have inflection points in m 2 and j3. In the 
region /3 > 2.2 the line separates the unconfined phase with vector mesons 
from the phase where confinement may occur. 
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In the region of weak coupling (large f3) and above the transition line 
we have evaluated the energy of the vacuum £ and the order parameter 
€ per site for m 2 = 0.1... 8. We have used a cubic lattice with size 8 
and 12 4 . No difference could be spotted in the results for the two choices. 
We have taken 10 3 measures separated by 15 updatings. For (3 = 10 we 
have performed 10 4 measures in order to reduce the statistical errors, but 
no appreciable improvement was observed. Some results are given in Figsj3] 
and |4l It is very interesting that both these global quantities converge for 
large m 2 to values independent from /3. From eq. (|30p one would expect 
a strong dependence from the mass. From both figures we see that near 
the transition line the dependence on f3 becomes strong. In Fig. 0]we show 
a simple fit of the data at f3 = 3, 10 with polynomials of second degree in 
ln(m 2 ). 



6 The Fit 



In the present Section we try a fit of the analytic results of Section [2] to the 
data obtained by simulation in Section [5j We use a simplified form of the 
mapping (|3|) 



9 = 

M 2 
A 2 
Ma 



4 

1 



m 

W) 



(34) 



The results of Section [2] are 

'(2) 3 



M 

(0) 
M 

(2) 
M 



_ 2 a 4 M 4 
2 9 jA^f 



8.625 



In m 2 - Ins - 4.939 



5.943 



+ e 



3 a 4 M 4 , 
3 , 



8.625 



0.348 In m 2 - In s - 2.954 + 0.232 



3 2 a 4 M 4 



(47T) 



8.625 



2 ( In m 2 - In s - 4.439 ) + 12.386 



(35) 



The parameters are fitted according to the following steps: i) first we fit the 
energy, since it is 0(g 2 ); ii) then we enter the parameters obtained by i) 
(aM and In s) into Cm and we perform the best fit for 



A 

aM 



-AI 



(2) 



M 



(36) 
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Figure 3: Energy per site at /3 = 3, 5, 6, 10, 20 as function of m 2 
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Figure 4: Order parameter £ per site at /3 = 3, 4, 5, 10, 20 as function of m 2 . 
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on a range of m 2 far from the transition line. A is a new free parameter in 
accordance to eq. (|34j> . 

The strategy is dictated by the fact that the direct fit of the order parameter 
does not determine properly the values of aM and Ins, since these last 
parameters have weak effect in comparison with the O(g ) term. Thus it is 
better to work at first with a 0(g 2 ) quantity: the vacuum energy. However 
the results of the first step depend strongly on the range of m 2 ; therefore 
we choose the interval that yields a common value for In s according to the 
parameterization of eq. (|34p . The interval is reported in table [1] in the 
second column. 



6.1 Fit of the Vacuum Energy 

We give few examples for (3 = 3, 4, 5, 6, 7, 8, 9, 10, 20. In Table [1] we list the 
results of the fit. It should be noticed that aM is almost constant with some 
small dependence on (3. x 2 is very small. Figs. and give a reasonable 
account of the match between lattice gauge theory and the field theory 
calculation with our subtraction procedure. 

6.2 Combined Fit of the Vacuum Energy and the Order Pa- 
rameter 

The second step allows us to fit the order parameter. Thus we find the value 
of A. The result of the procedure is listed in Table [2j In Figs. ([7]) and ([8]) 
we depict the results of the fit procedure for the values j3 = 3, 10. The plot 
include also the value of the Wilson action of the vacuum £ — (£. The figures 
show explicitly that the fit procedure intentionally excludes the region near 
the transition line (at m 2 ). 

The numerical results show that the field theory two-loop calculation of 
the energy of the vacuum and the order parameter provides a good fit of the 
lattice data on the surface given by the parametric equations 

2 4 
— — 

Ma ~ 2.3 

^-3m 2 , (37) 

which depart substantially from the naive identification m = Ma. This 
seems to be the unavoidable price for a mapping of the continuum theory 
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Figure 5: Fit of the field theory prediction (eq. (|35p ) to the lattice data for 
j3 = 3. The parameters are aM, Ins, eo- 
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Figure 6: See caption of Fig. [5] for /? = 10. 
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equipped with more scales than the lattice can accommodate. 

7 Conclusions 

We have investigated the deconfined phase of a massive Yang-Mills model. 
In particular we have considered the weak coupling limit (large f3 and m? S> 
m 2 ) of two global observables: the vacuum energy density and the order 
parameter density. The lattice simulation data (dependent on (3 and m 2 ) 
have been fitted with the analytic two-loop calculations of the continuum 
field theory (dependent on the mass M, the ultraviolet cut-off A and the 
lattice spacing a). The fit turns out to be very promising both in supporting 
the conjecture of the lattice as a regulator for nonlinear gauge theories and 
in the comprehension of the lattice parameters in phenomenology. 

The outlook is for a comparison of the lattice data on the energy gaps 
with the radiative correction of the masses (self-energies). Moreover the 
lattice simulations might provide quantitative predictions near the phase 
transition line |10j . 
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Figure 7: Fit of the field theory predictions (eqs. ([35]) . ([36]) ) to the lattice 
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Table 1: The fit of energy of the vacuum according to (|35|) . Ins is defined 
in eq. ([SI)- 





fit range 


aM 


In s 


x 2 


3 


[0.5:8.5] 


2.556 ± 0.010 


-1.009 ± 0.047 


4.3e-06 


4 


[0.218:8.5] 


2.475 ± 0.028 


-0.99 ± 0.16 


7.8e-05 


5 


[0.165:8.5] 


2.439 ± 0.039 


-1.07 ± 0.23 


l.le-04 


6 


[0.147:8.5] 


2.414 ± 0.043 


-1.13 ± 0.26 


9.3e-05 


7 


[0.125:8.5] 


2.416 ± 0.048 


-1.28 ± 0.28 


9.8e-05 


8 


[0.115:8.5] 


2.380 ± 0.047 


-1.17 ± 0.29 


7.2e-05 


9 


[0.109:8.5] 


2.383 ± 0.050 


-1.28 ± 0.30 


6.7e-05 


10 


[0.106:8.5] 


2.318 ± 0.046 


-0.97 ± 0.31 


4.0e-05 


20 


[0.07:8.5] 


2.313 ± 0.036 


-1.03 ± 0.25 


1.6e-05 



Table 2: The fit of the order parameter according to (|35p . A is defined in 
eq. ®. 



p 


fit range 


A 


3 


[1.0 


8.5] 


2.151 ± 0.011 


4 


[1.0 


8.5] 


2.2357 ± 0.0089 


5 


[1.0 


8.5] 


2.2774 ± 0.0081 


6 


[1.0 


8.5] 


2.2906 ± 0.0082 


7 


[1.0 


8.5] 


2.3067 ± 0.0079 


8 


[1.0 


8.5] 


2.3155 ± 0.0078 


9 


[1.0 


8.5] 


2.3237 ± 0.0077 


10 


[1.0 


8.5] 


2.3273 ± 0.0075 


20 


[1.0 


8.5] 


2.3492 ± 0.0051 
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